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and Athanasios V. Vasilakos 

Abstract 

In this paper, we propose an analog network coding scheme with differential modulation (ANC- 
DM) using amplify-and-forward protocol for bidirectional relay networks when neither the source 
nodes nor the relay knows the channel state information (CSI). The performance of the proposed 
ANC-DM scheme is analyzed and a simple asymptotic bit error rate (BER) expression is derived. 
The analytical results are verified through simulations. It is shown that the BER performance of 
the proposed differential scheme is about 3 dB away from that of the coherent detection scheme. 
To improve the system performance, the optimum power allocation between the sources and the 
relay is determined based on the simplified BER. Simulation results indicate that the proposed 
differential scheme with optimum power allocation yields 1-2 dB performance improvement over 
an equal power allocation scheme. 

Index Terms 

Differential modulation, bi-directional relaying, analog network coding, amplify-and-forward 
protocol 

I. INTRODUCTION 

Bi-directional relay communication has attracted considerable interest recently [OQ-O, 
and various bi-directional relay protocols for wireless systems have been proposed [DQ-[[3]|. 
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In El B3, the conventional network coding scheme is applied to the bi-directional relay 
network. Two source nodes transmit to the relay, separately. The relay decodes the received 
signals, performs binary network coding, and then broadcasts network coded symbols back 
to both source nodes. However, this scheme may cause irreducible error floor due to the 
detection errors which occur at the relay node. 

In fl3j — flSD , an amplify and forward based network coding scheme, referred to as the analog 
network coding, was proposed. In this scheme, both source nodes transmit at the same time so 
that the relay receives a superimposed signal. The relay then amplifies the received signal and 
broadcasts it to both source nodes. Analog network coding is particularly useful in wireless 
networks as the wireless channel acts as a natural implementation of network coding by 
summing the wireless signals over the air. 

Almost all existing works in bi-directional relay communications using analog network 
coding assume that the sources and the destination have perfect knowledge of channel state 
information (CSI) for all transmission links. As a result, coherent detection can be readily 
employed at either sources or relay, or both [OQ-|0. In some scenarios, e.g. the slow fading 
environment, the CSI is likely to be acquired by the use of pilot symbols. However, when 
the channel coefficients vary fast, channel estimation may become difficult. In addition, 
the channel estimation increases computational complexity in the relay node and reduces 
the data rate. Moreover, it would be difficult for the destination to acquire the source-to- 
relay channel perfectly through pilot signal forwarding without noise amplification. Hence, 
differential modulation without need of any CSI would be a practical solution. 

In a differential bi-directional relay network, each source receives a superposition of 
differentially encoded signals from the other source, and it has no knowledge of CSI of both 
channels. All these problems present a great challenge for designing differential modulation 
schemes in two-way relay channels. In [6], differential receivers for differential two-way 
relaying were presented using analog network coding. These non-coherent schemes were 
realized by averaging the channel coefficients of the substraction of adjacent received signals. 
However, the approaches may result in more than 3 dB performance loss compared to the 
coherent schemes due to instantaneous detection errors. 

In this paper, we propose an analog network coding scheme with differential modula- 
tion (ANC-DM) using amplify-and-forward protocol for bidirectional relay networks so that 
the CSI is not required at both sources and the relay. The performance of the proposed 
ANC-DM scheme is analyzed and a simple asymptotic bit error rate (BER) expression is 
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derived. The analytical results are verified through simulations. They show that the proposed 
differential scheme is about 3 dB away compared to the coherent detection scheme. To 
improve the system performance, the optimum power allocation between the sources and 
the relay is determined based on the provided simplified BER. Simulation results show that 
optimum power allocation yields 1-2 dB performance improvement over an equal power 
allocation scheme. 

Note that unlike [6], the destination realizes differential detection by subtracting away its 
own contribution in the received signals based on the power estimation. Besides, our scheme 
requires linear complexity, while the detectors in |[6l are much complex and non-linear which 
may calculate the modified Bessel function of second kind. Note also that the application of 
differential modulation for the bi-directional relaying in HI [0 with digital network coding 
is relatively straightforward as the received signals at the relay from the sources can be 
decoded separately due to the use of orthogonal transmissions. However, for amplify-and- 
forward relaying O-O, the relay receives a superposition of differentially encoded signals 
from the sources which makes the final detection at the source difficult. 

The rest of the paper is organized as follows: In Section II, we describe the proposed 
differential modulation scheme. Theoretical analysis is given in Section III. Section IV 
presents the optimal transmit power allocation between the sources and the relay. Simulation 
results are provided in Section V. In Section VI, we draw the main conclusions. 

Notation: Boldface lower-case letters denote vectors, (•)* stands for complex conjugate, 
(•) T represents transpose, E is used for expectation, Var represents variance, ||x|| 2 = x^x, 
and Jft(-) denotes real part. 

II. Differential Modulation for bidirectional relay networks 

A. Differential Encoding 

We consider a three-node bi-directional relay network consisting of two source nodes, 
denoted by Si and S 2 , and one relay node, denoted by R. All nodes are equipped with one 
antenna and operate in a half-duplex way so that the complete transmission can be divided into 
two phases, as shown in Fig. [TJ In the first phase, both source nodes simultaneously send the 
differentially encoded information to the relay, and in the second phase, the relay broadcasts 
the combined signals to both sources. Let c\(t)EA denote the symbol to be transmitted by 
the source S\ at the time t, where A represents a unity power M-PSK constellation set. In 
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the differential modulation bi-directional relay system, the signal Si(t) sent by the source Si 
is given by 

s 1 (t) = s 1 (t-l)c 1 (t), ci(t)eA (1) 
Similarly, the signal transmitted by S 2 at the time t is given by 

s 2 (t) = s 2 (t-l)c 2 (t), c 2 (t)eA (2) 

B. Differential Decoding 

In the bi-directional relayed transmission, the source nodes first broadcast the information 
to the relay. For simplicity, we assume that the fading coefficients are constant over one 
frame of length L, and change independently from one frame to another. The received signal 
in the relay at time t can be expressed as 

Vr{t) = y/PihiSi{t) + ^p~ 2 h 2 s 2 {t) + n r (t). (3) 

where pi and p 2 represent the transmit power at Si and £2 respectively, hi and h 2 are 
the Rayleigh fading coefficients with zero mean and unit variance between Si and R, and 
between S 2 and R, respectively, n r (t) denotes zero mean complex Gaussian random variable 
with two sided power spectral density of Nq/2 per dimension, and we furthermore assume 
Si, S 2 , and R have the same noise variance. 

In the second phase, the relay R amplifies y r (t) by a factor j3 and then broadcasts its 
conjugate, denoted by y*(t), to both Si and S 2 . The corresponding signal received by Si at 
time t, denoted by yi(t), can be written as 

Vi(t) = Py/frhiy;(t) + m(t) 

= fisl(t) + iys* 2 (t)+wi(t), (4) 

where (3 = (pi\hi\ 2 + p 2 \h 2 \ 2 + iVo)~ 5 , p r represents the transmit power by the relay, /j = 
l3y/pip r \hi\ 2 , v = (3^p 2 p r hih* 2 , andwi(t) = /3y/p^hin*(t)+rii(t). Note that unlike traditional 
ANC schemes O-O, y*(t) is transmitted from the relay, which obviously yields /x > 0. 
The reason of doing this is to make the receiver easily estimate fi and then subtract ixs\{t) 
in © for differential detection. 

As the relay has no CSI, the normalization factor (3 has to be obtained indirectly. We may 
rewrite the received signals in © in a vector format, given by 

y r = a/PiMi + VpHh 2 s 2 + n r , (5) 
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where y r = [y r (l), . . . , y r (L)] T , s x = [si(l), . . . , si(L)] T , s 2 = [s 2 (l), • • • , s 2 (L)] T , and 
n r = [n r (l), . . . ,n r (L)] T . To estimate the average receive power, we multiply the received 
signals by its Hermitian transpose as 

||y r || 2 = pi|/ii| 2 sfsi + p 2 \h 2 \ 2 s% s 2 + 2^f 2 y\{h\h 2 s*s 2 } 

+ 2^K{hxn^ Si} + 2 v /^H{/i 2 nf s 2 } + nf n r (6) 
By taking the expectation of ©, /3 can be then approximated at high SNR by 

where E{sf Si} = E{sf s 2 } = L, E{nf n r } = LJV , and E{sf s 2 } = E{nf si} = E{nf s 2 } = 
0. 

Similarly, the received signal at S 2 can be calculated as 

y 2 (t)=p^;h 2 y*(t)+n 2 (t), (8) 

As Si and S 2 are mathematically symmetrical, as shown in © and ©, for simplicity, we 
in the next only discuss the decoding as well as the corresponding theoretical analysis for 
signals received by Si. 

Recalling the differential encoding process in ©, © can be further written as 

yi(t) = fisl(t) + vs* 2 (t) + wi(t), 

= lM8l(t) + us* 2 (t - l)c£(t) + Wi{t), (9) 

Obviously, since si(t) is known, to decode c 2 (t) in ©, /xs*(£) should be removed. To do 
this, we have to first estimate where /i > 0. 

Similar to (flOl) . the received signal in © can be expressed in a vector form as follows 

y x = //si + vs 2 + wi, (10) 

where y x = [j/i(l), . . . , and Wi = [wi(l), . . . , Wi(L)] T . At high SNR, we may 

approximately obtain 

/ i 2 + |z/| 2 ^yfy 1 /L. (11) 

Since the source node Si can retrieve its own information si(t — 1) and Ci(t), based on (OQ) 
and ©, we have the following transformation 

yi(t)±cl(t)yi(t-l)-yi(t) 

= us* 2 (t - 1) (ci(t) - c 2 (t))* + wi(t), (12) 
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where w\(t) = cl(t)wi(t — 1) + Wx{t). Then, \u\ 2 can be approximately calculated in a 
similar way as (fTTI) 

~H~ 

\u\ 2 « illl (13) 

where y l = [yi(l), ...,yi(L- 1)] T , R[\s 2 (t - 1)| 2 ] = 1 , and the calculation of E[|ci(t) - 

c 2 (t)| 2 ] is given in Section-II-C. As ll is positive, by combining ©, ([111) and (PT3T) . we have 

f VA, A > 
/i « { (14) 
[0, A < 

where A = ^ili. — tw^tts — tttt^t- Note that in low SNR, A could be negative due 

L LE[|s 2 (t-l)| 2 |ci(t)-c 2 (t)| 2 ] b 

to the noise effect, and thus we set fj, ~ instead. The estimation method given in (fl4l) is 
evaluated in Fig. [2l 

By subtracting /is*(t), © can be further written as 

y'i(t) -Vi(t) -fia^t) 

= us* 2 (t-l)c* 2 (t)+w 1 (t) 

= (y[(t - 1) - Wl (t - 1)) c* 2 (t) + Wl (t), (15) 

Finally, the following linear decoder can be used to recover c 2 (t) 

c 2 (t) = arg max Re {y'My'*^ - l)c 2 (i)} . (16) 
c 2 (t)e.A 

And Ci(t) can be differentially decoded in a similar way by the source S 2 . Note that in 
comparison to traditional differential modulation 0, the extra complexity comes from the 
fx estimation in (PT4l) . which is linear and only comprises a few number of additions and 
multiplications. But the receiver in (6]| requires complicated computations, such as the zeroth- 
order modified Bessel function of the second kind. 

By ignoring the second order term, the corresponding SNR of the proposed differential 
detection scheme can be written as 
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where Var{ Wl (t)} = f5 2 p r N \h l \ 2 + N , ^ ^ Pl /N , ^ 2 = p 2 /N , and </v = p r /N { 
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Alternatively, if coherent detection is used, under the assumption of hi, h 2 and N available 
at Si, after subtracting pLSi(t)* from y 1 (t) in ©, the corresponding SNR can be calculated 
as 

A H 2 

7c Var{wx(t)} 
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| 2 + ^2i 




2 + l 



(18) 



By comparing (fTTT) and ([TBI , we can easily obtain 

7d ~ ^ (19) 

which clearly indicates that the differential detector in (PT6l) suffers around 3 dB performance 
loss compared to the coherent scheme. 

C. 77ze Calculation ofE[\ci(t) — c 2 (t)| 2 ] in ( Q 

From (fT3l . the average power of ci(t) — C2(i) needs to be calculated. When Af-PSK 
constellations are applied, the number of symbols produced in the new constellation by 
ci(t) — c 2 (t) is finite. Hence, it is easy to derive the average power of the new constellation 
sets. Note that the value of ci(t) — c 2 (t) can be equal to zero, which may affect the estimation 
accuracy in (fT3l . 

In order to overcome this problem, we may properly choose a rotation angle for the symbol 
modulated in source S 2 by c 2 (t)e~^ 6 ', ensuring that ci(t) — C2(t) in (fT3l is nonzero. For a 
M-PSK constellation, the effective rotation angle is in the interval [— 7r/M, ir/M] from the 
symmetry of symbols. For a regular and symmetrical constellation, the rotation angle may 
be simply set as 9 = ir/M. Similar approach may be used to generate the rotation angle for 
other types of constellations. 

Here, we give two examples on how to compute the average symbol power: 

1. Supposing BPSK constellation { — 1,1} is used, we have ci(t) — 02(f) E {—2,0,2}. 
Hence, the calculation of average power in the new set is straightforward. 

2. Supposing Si uses the BPSK set { — 1, 1}, by constellation rotation, S 2 can use {—j,j}- 
And we can get ci(t) — c 2 (t) E {— 1 — j, — 1 +j, 1 — j, 1 Hence, it is also easy to derive 
the average power in the new constellation set. 
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III. Performance Analysis 

For simplicity, we in this section analyze the BER performance using BPSK for the 
proposed ANC-DM scheme, and we assume p 1 = p 2 = p s and p s = Xp r , where A > 0, 
and thus, ipi = ip 2 = 4>s — AV+ (TTVT) can be rewritten as 
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where ip' r = (1 + \)ip r . 

Let X = 7^, and the BER for BPSK modulation can generally be expressed as 



BER = E [Q(2X)} 



exp(—x)Fx(x) 



x 



dx, 



(21) 



2^ Jo 

where Q(-) is the Gaussian-Q function, Fx(x) is the cumulative distribution function (CDF) 
of X. The right side of the equation can be readily obtained by integration by parts. The 
above expression is useful as it allows us to obtain the BER directly in terms of the CDF of 
X. 

By using a general result from [[8]|, the BER in (TfTI) can be approximated in the high SNR 
regime by considering a first order expansion of the CDF of X. Specifically, if the first order 
expansion of the CDF of X can be written in the form 



Fx(x) 



ax 



N+l 



x n+ \n + i) 



+ o(x 



,e > 0, 



(22) 



where A represents the average transmit SNR. At high SNR, the asymptotic BER is given 
by GO 



BER 



aT(N + 



2^X N+1 (N+1) 



+ o((A)"^ +1 )). 



(23) 



The PDF of X can be obtained with the help of [|9]| 

8(1 + \fx exp (-2(1 + X)^- 1 + ip'-^x) 



Px(x 



1ps + 1p' r 



4(1 + X)x 



+2K 



4(1 + X)x 



U{2{l + X)x) 



(24) 

where K (-) and K\{-) are the zeroth-order and first-order modified Bessel functions of 
the second kind, respectively, and U(-) is the unit step function. Note that the exact BER, 
which is complicated in computation, does not have a closed-form solution. However, at high 
SNR, when z approaches zeros, the K\(z) function converges to 1/z ifTOl . and the value of 
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the K (z) function is comparatively small, which could be ignored for asymptotic analysis. 

Hence, Px(x) in (1241) can be approximated as 

p , > _ 2(1 + A)(-fc + exp (-2(1 + A)(fe' + 

' ( ' 

For the ANC-DM, the CDF of destination SNR ^ d can be approximated as 

F x (x) « 1 - exp (-2(1 + A)(C' + VO) 

« 2(1 + A)^; 1 + O + o(x 1+£ ). (26) 

Finally, comparing (g6]) with and (1231) . the asymptotic BER of ANC-DM at high SNR 
can be approximated as 

= (i + X)(^ + ^) ^ 



IV. Transmit Power Allocation 

In this section, we discuss how to allocate power to both sources and the relay subject to 
total transmission power constraint. It can be seen from (1271) that the asymptotic BER of the 
proposed differential modulation scheme depends non-linearly upon p s and p r . Hence, when 
the total transmit power is fixed, 2p s + p r = p, the power allocation problem over Rayleigh 
channels can be formulated to minimize the asymptotic BER at high SNR in (|27T) 

min BER 

S.t. 2p s +p r =P (0<p s <p, 0<p r <p), (28) 

where we assume p± = P2 = p s and p s = Xp r . 

The power allocation problem is to find p s such that the BER in (1271) is minimized subject 
to the power constraint by solving the following optimization problem 

C( Ps ) = BER + ^2p s +p r -p), (29) 

where £ is a positive Lagrange multiplier. The necessary condition for the optimality is found 
by setting the derivatives of the Lagrangian in (1291 with respect to p s and p r equal to zero, 

respectively. Reusing the power constraint, we can calculate at high SNR that 

p 

Pr = |, (30) 
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which indicates the power allocated in the relay should be equal to the total transmit power 
at both sources in order to compensate the energy used to broadcast combined information 
in one time slot. 

V. Simulation Results 

In this section, we provide simulation results for the proposed ANC-DM scheme. We 
also include corresponding coherent detection results for comparison. All simulations are 
performed for a BPSK modulation over the Rayleigh fading channels. The frame length is 
L = 100. For simplicity, we assume that 2p s + p r = p = 3, and Si, S2 and R have the same 
noise variance N . The SNR ip s can be then calculated as ip s = p s /N . 

Fig. |2] shows the simulated BER performance for differential and coherent ANC schemes in 
bi-directional relaying without using constellation rotation. Equal transmit power allocation 
is applied: p s = p r = p/3. It can be seen that the differential scheme suffers around 3 
dB performance loss compared to the coherent ANC scheme, which has been validated by 
(fl9l ). We also include the Genie-aided result by assuming that ji is perfectly known by the 
source such that traditional differential decoding can be performed. It shows from the results 
that there is almost no performance loss using the estimation method in (fl4l) which clearly 
justifies the robustness of the proposed differential decoder. 

It is worthwhile mentioning that, in ((J), it uses similar Genie-aided result as a benchmark 
as well. The major difference is that the detectors in [6] have much inferior performance than 
the Genie-aided result, and it has about 6dB performance loss in comparison to the coherent 
detection scheme. However, our proposed detection algorithm has comparable performance 
with the genie-aided result, and only has 3dB performance loss than the coherent detection 
results. This clearly indicates that our proposed method outperforms the differential detectors 
in (6]|. The main performance loss in [[61 is due to that uncoherent detection approach is 
employed by statistically averaging off the impact of channel fading coefficients ignoring 
the instantaneous channel state information. But, our method utilizes differential detection 
relying on the operation with the previous received signals which could be more adaptive to 
variation of the channels. 

In Fig. [31 we compare the analytical and simulated BER performance of the proposed 
differential modulation scheme. Equal transmit power allocation is also applied by setting 
Ps = Pr — p/3- From the figure, it can be observed that at high SNR, the analytical BER 
derived by (1271) converged to the simulated result, which justifies the validation of (1271) . 
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In Fig. HI we examine the BER performance of the proposed differential modulation 
protocol with power allocation by setting p s — p/A and p r — p/2 subject to the total power 
constraint. From Fig. |H it can be observed that with optimal power allocation, the proposed 
scheme obtains about 2 dB performance gain in comparison with the equal power allocation 
scheme at high SNR. We also compare the result by another power setting: p s = OAp and 
p r = 0.2p, and inferior result can be again observed compared to the optimal power allocation. 

In Fig. [5l we plot the BER curves in terms of A = p s /p r defined in (l20l using different 
noise variance N with the asymptotic BER constraint in (1271) . With the power constraint 
2p s + Pr = P, the SNR can be therefore calculated by ip s = ( 2 a+iw > wnere p = 3. It shows 
that best performance is obtained when A = 0.5. In other words, p s = p/A and p r = p/2 is 
the optimal power setting between the sources and the relay, which further verify the power 
allocation strategy in (l30l) . From the figure, we can also see that the asymptotic BER is very 
close to the simulated results for various SNR values, and they result in the same power 
allocation solution of p s = p/A and p r = p/2. 

In Fig [61 we examine the BER results of the proposed differential modulation scheme 
without and with using constellation rotation, where the signal constellation used by S\ is 
rotated by n/2 relative to that by 5*2. It can be observed that the new result has very similar 
with the curve without rotating constellations. This indicates that using constellation rotation 
may not give system any gains given large frame length. 

Fig. [7J shows the average normalized mean square error (MSE) of /i estimation in (PT41) as a 
function of the SNR, where the average normalized MSE is calculated as - ^ k= ^^ 
and Jl(k) is the estimate of p.(k). It can be observed that the estimation is quite accurate, 
particularly at high SNR. 

VI. Conclusions 

In this paper, we have proposed a simple differential modulation scheme for bi-directional 
relay communications using analog network coding when neither sources nor the relay has 
access to channel state information. Simulation results indicate that there exist about 3 dB 
loss compared to the coherent detection scheme. Analytical BER is derived to validate the 
proposed method. In addition, based on the asymptotic BER at high SNR, an optimal power 
allocation between the sources and the relay was derived to enhance the system performance. 
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Fig. 1. Block diagram of the proposed ANC-DM scheme. 




Simulated BER performance by differential and coherent detections, where p s — p r = p/3. 




Fig. 4. Simulated BER performance by the proposed differential scheme using transmit power allocation. 
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Fig. 5. Simulated BER performance by the proposed differential scheme using transmit power allocation in term of 
A = Ps/pr, with different noise variance No and ip 3 — 



(2\+l)N ■ 
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Fig. 6. Simulated BER performance by the proposed differential scheme with and without using constellation rotation, 
where p a — p r = p/3. 




Fig. 7. The average normalized mean square error (MSE) of fi estimation in d 14b as a function of the SNR. 



